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The thermodynamic properties of the bilayer Sr3Ru2C>7 at very low temperatures are investigated 
by a realistic tight-binding model with the on-site interactions treated at the mean-field level. Due 
to the strong spin-orbit coupling, the band structure undergoes a significant change in Fermi surface 
topology as the external magnetic field is applied, invalidating the rigid band picture in which the 
Zeeman energy only causes chemical potential shifts. In addition, since Sr3Ru207 is a t2 g active 
system with unquenched orbital moments, the orbital Zeeman energy is not negligible and plays an 
important role in the phase diagram on the magnetic field orientation. We find that both the total 
density of states at the Fermi energy and the entropy exhibit a sudden increase near the critical 
magnetic field for the nematic phase, echoing the experimental findings. Our results suggest that 
extra cares are necessary to isolate the contributions due to the quantum criticality from the band 
structure singularity in this particular material. The effects of quantum critical fluctuations are 
briefly discussed. 

PACS numbers: 72.80.Ga,73.20.-r,71.10.Fd 



I. INTRODUCTION 

Recently, a great deal of attention has been paid to 
the bilayer ruthenate compound Sr3Ru20y with various 
interesting physical properties. It was first considered 
as a field-tuned quantum critical state having metamag- 
netic transitions around 8 Teslai^— . Later, in the ultra- 
pure single crystal it has been found that the meta- 
magnetic quantum critical point is intervened by the 
emergence of an unconventional anisotropic (nematic) 
electronic stated, stimulating considerable theoretical 
efforts^—. Sr3Ru2 07 is a metallic itinerant system with 
the active ^-orbitals of the Ru sites in the bilayer RuC>2 
(afe) planes. At very low temperatures (~ IK), it starts as 
a paramagnet at small magnetic fields. Further increas- 
ing field strength leads to two consecutive metamagnetic 
transitions at 7.8 and 8.1 Tesla if the field is perpen- 
dicular to the a6-plane. The nematic phase is observed 
between these two transitions, identified by the observa- 
tion of anisotropic resistivity without noticeable lattice 
distortions. 

This nematic phase in Sr3Ru2 07 can be understood as 
the particle-hole channel Fermi surface instability of the 
Pomeranchuk-typc 3 . It is a mixture of the (i-wave Pomer- 
anchuk instabilities in both density and spin channels^ 9 -, 
though the microscopic origin of these instabilities re- 
main controversial. Different microscopic theories have 
been proposed based on the quasi-lD bands of d xz and 
d ys M^^., and based on the 2d-band of d x ^^L. In 
ou r 14 ' 18 and Raghu et. a/.'s 15 theories, the unconven- 
tional (nematic) magnetic ordering was interpreted as or- 
bital ordering among the d xz and dj^-orbitals. In partic- 
ular, in Ref. [18| a realistic tight-binding model with the 
important features including the i2g-orbitals, the bilayer 
splitting, spin-orbit coupling, and the staggered rotations 
of RuO octahcdra, has been derived and shown to repro- 
duce accurately the results of the angle-resolved photon 
emission spectroscopy (ARPES) 20 - and the quasiparticle 



interference in the spectroscopic imaging scanning tun- 
neling microscopy (STM) 21 -. 

On the other hand, orbital physics is an important sub- 
ject in transition metal oxides with active ri-orbital a 22 ' 23 . 
Recently the research of orbital physics has been fur- 
ther stimulated due to the observations of anisotropic 
properties in iron-pnictide a 24 ' 25 which are consistent with 
orbital ordering 2 ^— . Generally speaking, a transition 
metal oxide with active i2g-orbitals could exhibit or- 
bital ordering in the quasi-lD bands under certain con- 
ditions. In the particular case of Sr3Ru20y, it has been 
argue d 14 ' 15 that at zero magnetic field, although the sys- 
tem does not show any order, it is near the instabilities 
of ferromagnetic and orbital (nematic) orderings. Ap- 
plying the magnetic field pushes the system closer to the 
van Hove singularities which largely raises up the density 
of states at the Fermi energy T>(ep) and results in both 
the metamagnetic and nematic transitions. It should be 
noted that the Fermi liquid states always dominate the 
low temperature physics of this system because the re- 
sistivity was found to follow T 2 law in both isotropic 
and nematic phases. This suggests that quantum criti- 
cality is not expected to play much role in both isotropic 
and nematic phases at low temperatures. Nevertheless, 
quantum critical behavior appears at finite temperatures 
above the nematic phase, which exhibits non-Fermi liquid 
properties. 

The study of quantum critical states has attracted 
intensive research activities over past few decade s 29 ' 30 . 
Materials possessing quantum critical points are usually 
triggered by tunable parameters, for example, external 
pressure and magnetic field. When approaching quan- 
tum critical points, not only the zero temperature elec- 
tronic ground state but also the low-temperature proper- 
ties are modified radically. Pioneering theoretical works 
done by Herta^i and extended by Millis 32 have addressed 
the physics near the ferromagnetic and antiferromagnetic 
quantum critical points, showing that the critical flue- 
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tuations are scale invariant in the vicinity of quantum 
critical points due to the divergence of the correlation 
length right at these points. This aspect enables the es- 
tablishment of universal classes based on the scaling the- 
ory. Experimentally, quantum critical states usually ex- 
hibit divergences in thermodynamic properties like spe- 
cific heat, entropy, etc. The critical exponents associated 
with these divergences are believed to be universal val- 
ues in the same universal class regardless the microscopic 
details. 

The influence of quantum critical fluctuations in 
Sr 3 Ru207 seems to be novel as well. Rost et. a?j 33 i 34 
measured the entropy and specific heat in ultra-pure sam- 
ples and found divergences near the metamagnetic transi- 
tions in both quantities. Although it is a common feature 
in a quantum critical state that the specific heat C di- 
verges as C/T = A[(H - H c ) / H c ]- a + B due to quantum 
fluctuations, the exponent of a is fitted to be 1 instead 
of 1 /3 as predicted by the Hertz-Millis theory 3 ^. The 
total density of states (DOS) measured by Iwaya et al£5. 
using the STM showed that the DOS at the Fermi en- 
ergy {T>(£p)) increases significantly under the application 
of the magnetic field, but the DOS at higher and lower 
energy does not change accordingly. This indicates that 
the Zeeman energy in this case does not simply cause 
a chemical potential shift so that the conventional rigid 
band picture can not explain this result. These findings 
have posted a challenge to understand the critical behav- 
ior in this material. 

In this article, we show that the tight-binding model 
derived by Arovas and two of us^ together with the on- 
site interactions treated simply at the mean-field level 
already gives rise to divergences in these properties with- 
out involving quantum fluctuations. Realistic features 
like multi-orbital bands, the bilayer splitting, the rota- 
tions of RuO octahedra, and spin-orbit coupling make 
this material very sensitive to small energy scales. Be- 
cause parts of the Fermi surface are close to the van Hove 
singularities, Fermi surface reconstructions in the exter- 
nal magnetic fields lead to a singular behavior in U(ef). 
This results in the divergences observed in the exper- 
iments mentioned above. Because of the strong spin- 
orbit coupling and the unquenched orbital moments, the 
Zeeman energy tends to reconstruct the Fermi surfaces 
rather than just provides a chemical potential shift. Our 
results suggest that the influence of quantum critical fluc- 
tuations will be masked if T>{cf) of the system exhibits 
a non-monotonic behavior with the tuning parameters, 
implying that a more careful analysis is required in order 
to distinguish the role of the quantum criticality in the 
bilayer Sr3Ru207. 

This paper is organized as follows. We will summarize 
the tight-binding model derived in Ref. [HI in Sect. Ql] 
The formalism of the mean-field theory will be presented 
in Sect. IIIIl The results for the cases of the magnetic 
fields parallel to and tilted away from the c-axis will be 
discussed in Sect. |IV] and [V] respectively. Conclusions 
will be given in Sect. IVI1 



II. THE TIGHT-BINDING MODEL 

The detailed band structure of Sr3Ru2 07 is compli- 
cated by the i2g-orbital structure (e.g. d xz , d yz , d xy ), the 
bilayer splitting, the staggered distortion of the RuO oc- 
tahedra, and spin-orbit coupling. We have constructed 
a detailed tight-binding Hamiltonian which gives rise to 
band structures in agreement with the ARPES data in 
a previous work 18 . We found that a difference of the 
on-site potential between the two adjacent RuO layers, 
Vbias, should be added-^ in order to fit the shape of the 
Fermi surfaces observed in the ARPES experiments^ . 
This term appears because ARPES is a surface probe 
and this bilayer symmetry breaking effect is important 
near the surface. Since we focus on the thermodynamic 
properties which are all bulk properties, Vbias is set to be 
zero in this paper. Below we will start from this model 
and refer readers to Ref. [l8| for more detailed informa- 
tion. 

The tight-binding band Hamiltonian Hq can be re- 
duced to block forms classified by k z = 0,ir correspond- 
ing to bonding and anti-bonding bands with respect to 
layers as: 

Ho = ho{k z = 0) + h (k z =n), (1) 
with ho(k z ) defined as 

h (k z ) = Y'$t ( h osfk z ) g\(k,k z ) \ 
K» A. k ,k m ,.\ g^ kz) h 0s (k + Q,k z )J k - k *> s 

k 

(2) 

where the spinor $t operator is defined as 

k,k z ,s 
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k+Q,s,k z ' k+Q,s,k z ' k+Q,-s,kJ , y ' 

k (k) annihilates an electron with orbital a and spin 

polarization s at momentum (k,k z ); Q = (tt,tt) is the 
nesting wavevector corresponding to unit cell doubling 
induced by the rotations of RuO octahedra; means 
that only half of the Brillouin zone is summed. Please 
note the opposite spin configurations s and -s for the 
d X z,d yz and d^-orbitals in Eq. [3J which is convenient 
for adding spin-orbit coupling later. 

The diagonal matrix kernels ho s in Eq. [5] are defined 

as 



h 0s {k,k z ) = A s (k) + Bi cos k z - fxl, (4) 
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where tj_ is the longitudinal inter-layer hopping for the 
d xz and d yz orbitals. A is the spin-orbit coupling strength 
which comes from the on-site spin-orbit coupling term 
as H so — ^^2 i Li -Sf, fi is the chemical potential; the 
dispersions for the d yz , d xz , and d xy bands in Eq. [5] are 
defined as 



uz 

XZ 
k 

XV 



-2t2 cos k x — 2t\ cos k y , 
-2ti cos k x — 2i 2 cos k yi 
-2t% ( cos k x + cos ky) — 
-2*5 ( cos 2k x + cos 2ky) 
-Atg sin k x sin k y , 



4*4 cos k x cos k y 

V X v 



(7) 



which includes longitudinal (t\) and transverse (£2) hop- 
ping for the the d xz and d yz orbitals, respectively, as 
well as are nearest neighbor (£3), next-nearest neighbor 
(£4), and next-next- nearest neighbor (ts) hopping for the 
d xy orbital. Following the previous LDA calculations 36 , 
V xy is introduced to account for the splitting of the d yz 
and d xz states relative to the d xy states. While symme- 
try forbids nearest-neighbor hopping between different 
tig orbitals in a perfect square lattice without the rota- 
tion of Ru octahedra, a term describing hopping between 
d xz and d yz orbitals on next-nearest neighbor sites (te) 
is allowed and put into the tight-binding model. 

The off-diagonal matrix kernel g(k, k z ) in Eq. [2] reads 



g(k,k z ) = G{k) -2B 2 cos k z , 



(8) 



where 



where a is the layer index, B = \B\, and the matrices 
L XjZ can be found in Ref. [1 

For 9^0, the extra Zeeman terms from x-component 
of L and S couple <£>t fe and $t fe . Defining <f)t ^ = 



matrix form as: 



I , the Zeeman term can be written in the 



Hzeernan = MB-B^ ^ (f)t fr H Z (0)(f)~ k , 
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(13) 
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with s = ±1 and 



H%(6) = sinfl x 




(15) 



III. THE MEAN-FIELD THEORY 
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and 



G{k) = 





2t mT j(k) 




-2t mT j{k) 





(9) 
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tot and t^ T describe the 



with j(k) = cosk x + cosk y . t 
intra- and inter-layer hopping between d xz and d yz in- 
duced by the rotations of RuO octahedra, providing the 
coupling between k and k + Q. 

When describing the Zeeman energy, we can choose 
the magnetic field B to lie on the xz plane and define 
9 as the angle between B and the c-axis of the sample 
without loss of the generality. Consequently, the Zeeman 
term becomes: 



TTOrbital 
Zeeman 



jjspin 

Zeeman 



Tjorbital , rrspin 
n Zeeman ' Zeeman' 

-fXB B ^ {Lz.ia cos 9 + L x . la sin 9) , 

i,a 

-2n B B {S a zla cos 9 + S a xla sin 9) , 



(11) 



In Sect. UH we have introduced the rather complicated 
band structure of Sr3Ru207. In this section, we intro- 
duce the minimal multi-orbital band Hubbard model for 
the nematic ordering and outline the mean-field proce- 
dure. 

In realistic band structures measured by ARPES 20 , 
there exists an additional <5-band arising from the d x 2_ y 2- 
orbital which is not covered by the current model. The 
particle filling in the ^-orbitals is not fixed. For the 
convenience of calculation, we fix the chemical potential 
\i = 0.94ti instead of fixing particle filling in the ti g - 
orbitals while changing magnetic fields and orientations 
in Sect. IIVI and Sect. [V] The corresponding fillings in- 
side the i2g-orbitals per Ru atom varies within the range 
between 4.05 and 4.06 in Figs. [21 [H EJ This treatment 
does not change any essential qualitative physics. 

The Hubbard model contains the on-site intra and in- 
ter orbital interactions as 



Hi, 



U 



E 



iaf ll ia\ 



I E • 



(16) 



where the Greek index a refers to the orbitals xz, yz and 
xy\ the Latin index a refers to the upper and lower layers. 
Throughout this paper, the parameter values are taken 
as U/ti — V/ti — 3.6. The other two possible terms 
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where 




FIG. 1: (Color online) The Fermi surfaces using the bilayer 
tight-binding model Eq. [2]with the parameters in unit of ti as: 
t 2 = 0.1ti,t 3 = ti,U = 0.2ti,t s = -0.06ti,t 6 = 0.1h,t± = 
0.6ti,ti NT — i INT = 0.1ti,X = 0.2ti,V xy — 0.3ti, and fi — 
0.94fi. The thick dashed lines mark the boundary of half 
Brillouin zone due to the unit cell doubling induced by the 
rotation of RuO octahedra. The Fermi surfaces of the bonding 
(k z = 0, black solid lines) and the anti-bonding bands (k z = 
it, red dashed lines) could cross since k z is a good quantum 
number. 

in the multi-band Hubbard interaction are the Hund's 
rule coupling and pairing hopping terms, which do not 
change the qualitative physics and are neglected. We 
assume the external B-field lying in the irz-plane with 
an angle 6 tilted from the z-axis. The occupation and 
spin in each orbital and layer are defined as follows: 

< = £«lM< a (z)>, 

s 

s 

s% a = I^KlW^W). (i7) 

No any other inter-layer interaction is considered and the 
non-interacting Hamiltonian has symmetric layers, thus 
the order parameters are the same for both layers, which 
have been numerically confirmed. From now on it will be 
assumed: 

< = n a , S« za =S« z . (18) 

With this property, k z remains a good quantum number 
in the resulting mean-field Hamiltonian. 

The standard mean-field decomposition of Hi nt leads 

to 

i,a,a s 

(19) 



The interaction parameters (U, V) are chosen such that 
no spontaneous magnetization occurs in the absence of 
the external magnetic field. Moreover, since the order pa- 
rameters are non-zero even without magnetic field, 
we require that the renormalized Fermi surface at zero 
field to be the one given in Fig. [T] As a result, in addi- 
tion to the optimized parameters obtained in our previ- 
ous worki£, we need to subtract the following term from 
Eq. [H 

H shift = EE^(°)^iw<a«, (2i) 

i,a,a s 

where 

W»(0) = i[/n Q (0) + V n/3 (°)' ( 22 ) 

and n a (0) is the occupation number in orbital a corre- 
sponding to the Fermi surfaces shown in Fig. [T] This is 
an effect of the renormalization of the chemical potential 
fi and V xy due to interactions. After putting all the pieces 
together, we finally arrive at the mean-field Hamiltonian 
as: 

H MF — Ho + H^f — Hghift + Hzeeman 
k k * 

The order parameters are computed self-consistently. 
It has been pointed ou t 14 ' 15 that the nematic phase can 
be identified as the orbital ordering between the d xz and 
<i az -orbitals, thus the nematic order parameter Af is de- 
fined as 

N = 2(n yz - n xz ). (24) 
The magnetization order parameter is defined as: 

M = 2^S a . (25) 

a 

The factor of 2 in Eq. [25] and [24] accounts for the double 
layers. 

IV. THE CASE OF THE PERPENDICULAR 
MAGNETIC FIELD (9 = 0) 

In this section the thermodynamic properties at low 
temperatures are calculated within the mean-field the- 
ory for the case of B \\ c, i.e., = 0. Most important 
features of these properties, including the evolution of 
the nematic ordering, the tunneling DOS, entropy land- 
scapes, and the finite temperature phase diagram under 
the application of the magnetic field B, can be reasonably 
reproduced and understood by the singular behavior of 
T>{ep) under the magnetic field. 
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FIG. 2: (Color online) The order parameters as a function 
of hbB for 8 = 0. The nematic phase (in the green area) is 
bounded by two magnetization jumps which corresponds to 
metamagnetic transitions. 



A. Nematic ordering and metamagnetic transitions 

The order parameters \A4 \ and J\f as functions of (jlbB 
for this case are shown in Fig. [2 There are three rapid 
increases in the magnetization, consistent with the ex- 
periment measurements of the real part of the very low 
frequency AC magnetic susceptibility at 7.5T, 7.8T and 
8. IT, respectively^. Experimentally, only the last two ex- 
hibit dissipative peaks in the imaginary part of the AC 
susceptibility, which characterize the first order meta- 
magnetic transitions. The first jump measured in the 
experiment is considered as a crossover. 

The nematic ordering develops in the area bounded by 
the last two magnetization jumps, reproducing the well- 
known phase diagram of the S^Ru^O?*. In particular, 
if we adopt the results from LDA calculation a 37 ' 38 that 
ti » 300 meV, we find that three jumps in the magne- 
tization appear at B rs 0.0032ti//is, 0.0053ii//is, and 
0.0059ii/£t B ~ 15.7 T, 26 T, and 29 T, which is within 
the same order to the experimental values of 7.5T, 7.8T, 
and 8. IT. This is an improvement compared to the results 
in previous theory calculations ^ 11 ' 15 ' 16 , in which the ne- 
matic ordering develops at much higher field strength 
fisB/h « 0.02 

The sensitivity to the small energy scale like the Zee- 
man energy is because a part of the Fermi surfaces, 
mostly composed of quasi-lD bands as shown in the 
yellow areas in Fig. |3] (b), is getting close to the van 
Hove singularities at (tt, 0) and (0, 7r). The evolution 
of the Fermi surface structures as increasing the B-ficld 
across the nematic phase boundaries is presented in Fig. 
[3] a (before), b (inside), and c (after) at nsB/ii = 
0.0048,0.00544,0.006, respectively. Before and after the 
nematic phases, the Fermi surfaces have the 4-fold ro- 
tational symmetry as exhibited in Fig. |3] (a) and (c). 
When the system is in the nematic phase, the Fermi 
surfaces only have 2-fold symmetry as expected. Par- 
ticularly, the nematic distortion is most prominent near 
(±7r, 0) and (0, ±7r) whose Fermi surfaces are dominated 
by the quasi-lD bands as shown in Fig. Efb), support- 
ing the mechanism of orbital ordering in quasi-lD bands 



driven by the van Hove singularities. 

It is worthy of mentioning that the onsite spin-orbit 
coupling H so — A Li ■ Si has important effects on 
the Fermi surface evolutions. H so hybridizes the oppo- 
site spins between quasi-lD bands d yZtXZ and the 2-D 
band d xy . As the spin Zeeman energy is present, the 
spin majority (minority) bands of d yz%xz couples to the 
spin minority (majority) band of d xy . Moreover, the 
orbital Zeeman energy provides more hybridizations be- 
tween quasi-lD d XZtVZ bands. Combined with the above 
two effects, the addition of the spin and orbit Zeeman en- 
ergies causes reconstruction of the Fermi surfaces rather 
than just chemical potential shifts. 

These results show that the complexity and the sensi- 
tivity of the Sr3Ru207 band structure can be captured 
very well by our tight-binding model with a reasonable 
quantitative accuracy. In the following, the same model 
will be used to further investigate some novel physical 
properties observed in experiments. 



B. Entropy landscape under the magnetic field 

One of the intriguing puzzles observed in Sr3Ru2 07 is 
the critical exponent of the divergence in entropy (as well 
as specific heat) when approaching the nematic region, or 
the quantum critical point. Although it is expected that 
specific heat C diverges as C/T = A[{H - H c )/H c ]~ a + 
B near the quantum critical point, the exponent a has 
been found to be l 33 ' 34 instead of 1/3 as predicted by the 
Hertz-Millis theory. However, it is noticed that T>(ep) is 
taken to be a constant in the Hertz-Millis theory while 
this material has a complicated Fermi surface evolutions 
under the magnetic field. Moreover, as mentioned in the 
introduction, the system is always in the Fermi liquid 
states at low temperatures inferred from the temperature 
dependences of the resistivity. Motivated by these two 
facts, it is then worthy of studying first the contribution 
from the band structure to the entropy before considering 
the quantum fluctuations. 

The entropy per Ru atom at the mean-field level can 
be evaluated by: 

S(B) = ^Ek'Ej In /(£,(£)) 

+(1 -/(£, (£)) ln(l -f(Ej (£)))], (26) 

where (3 — 1/ksT, f(Ej(k)) the Fermi-Dirac distribu- 
tion function, and Ej(k) is the jth eigenvalue of H MF (k) 
given in Eq. 1231 Because only the change of the entropy 
as a function of magnetic field near the nematic region 
is experimentally relevant, we plot in Fig. |4] the quan- 
tity: S(B)/T in unit of k 2 B jt\ at a low temperature of 
l/(/3ti) = 0.002 within the range of 0.0045 < n B B/h < 
0.007. Near the nematic region, S(B)/T increases first, 
being cut off inside the nematic region, and then de- 
creases after, which is consistent with the experiment 33 . 
Since these results are obtained at the mean-field level, 
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FIG. 3: Fermi surface evolution for the case of magnetic field parallel to c axis (a) before (fisB — 0.0048ti) (b) inside 
({ibB = 0.00544ii), and (c)after ([ibB = 0.006£i) the nematic phase. Significant changes in the Fermi surface topology under 
the magnetic field can be seen (see Fig. [T]for the Fermi surfaces at zero field). The nematic distortion is most obvious in the 
Fermi surfaces near (±7r,0) and (0, ±tt) as indicated by the yellow areas in (b). These parts of Fermi surfaces are composed 
mostly of quasi-lD bands, supporting the intimacy of nematic phase to the orbital ordering. 



this diverging behavior near the nematic region of the 
entropy landscape results from the singularity of T)(ep). 
As apparent in Eq. 1261 at very low temperature only the 
states near the Fermi surfaces give a sizable contribution 
to the entropy. Since the nematic region is driven by 
the sudden increase of V(cf), the entropy should also be 
enhanced near the nematic region. While it is generally 
expected that the quantum fluctuations near the critical 
point could contribute more entropies, our results demon- 
strate that entropy contributed from the mean-field level 
already gives diverging behavior due to the singularity of 
T>(cf) under the evolution of the magnetic field. 

To summarize, it is found that the singular behavior 
of T>(ep) already produces diverging behavior in both en- 
tropy and specific heat under magnetic field at constant 
temperature, although the critical exponent a can not be 
extracted from the current theory. Similar argument has 
been proposed in previous study 3 -^, in which the effect of 
a rigid band shift away from van Hove singularities in a 
perfect 1-d band is discussed. The quantum fluctuations 
seem to play a minor role in this case but become im- 
portant in the vicinity of the quantum critical point and 
at higher temperature. It is noted that the rigid band 
picture does not work neither when the system is doped 
by the substitution of La 3+ into the Sr 2+ site, which has 
been argued to result from the strong correlation effect 39 . 
We will leave these discussions in the conclusion section. 



C. Total density of states 

Iwaya et al. has measured the STM tunneling differen- 
tial conductance Mr in the S-field for Sr3Ru207, which 
corresponds to measurement of the DOS. It has been ob- 
served that while DOS at higher and lower energy does 
not change, the DOS at the Fermi energy (T>(cf)) in- 




0.005 0.006 0.007 

M B B/t l 

FIG. 4: (Color online) (a) The entropy landscape represetned 
by the quantity S(B)/T in unit of k%/ti within the range 
of 0.0045 < fiBB/ti < 0.007. A sudden increase near the 
nematic region (green area) is clearly seen. 

crease significantly under the application of the magnetic 
field, demonstrating the violation of the rigid band pic- 
ture. 

In our model the total DOS can be evaluated using: 
Ptot{u) = ^TrIm[G MF (fc, W )] 

k 

G MF (k,u) = (uj + if 1 -H MF (k)y 1 (27) 

where H MF (k) is given in Eq. [23]with the self-consistent 
order parameters and N is the total number of sites in 
the bilayer square lattices. 

The profiles of the total DOS at several different mag- 
netic field strength are plotted in Fig. [SJa), and clearly a 
rigid band picture does not apply here. T>(cf) (ptot{w — 
0) in the plot) increases significantly as the nematic phase 
is approached. This feature can also be directly un- 
derstood by the picture of Fermi surface reconstruction, 
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FIG. 5: (Color online) (a) The total DOS ptot(w) as a function 
of \xbB. The broadening factor r\ is set to be 77 = 0.002£i. 
ptot(ui) does not follow the rigid band picture and p«ot(w = 0) 
has a sudden increase near the nematic region. Inset: the total 
DOS at zero field for a wider range of \w\/t\ < 0.4. The peaks 
corresponding to van Hove singularities are near ui/t\ ~ —0.2 
and uj/t\ m 0.35. (b) The DOS of the quasi-lD bands. Inset: 
the DOS of quasi- ID bands at zero field for a wider range 
of |w|/ti < 0.4. Only the peak around u/ti w 0.35 remains, 
meaning that this peak is due to the van Hove singularities 
in quasi- ID bands. 



gion is passed. To resolve this discrepancy, several realis- 
tic features need to be considered before comparing our 
calculations with the STM results. Since the STM is a 
surface probe and the surface of the material is usually 
cleaved to have the oxygen atoms in the outermost layer, 
there is an oxygen atom lying above each uppermost Ru 
atom. Consequently, the tunneling matrix element will 
be mostly determined by the wavefunction overlaps be- 
tween the p-orbitals of the oxygen atom and the d-orbitals 
of the Ru atom, resulting in a much smaller matrix ele- 
ment for d xy orbital compared to d xz ^ yz orbitals 18 . The 
minimal model to take this effect into account is to ex- 
tract the DOS only from the quasi-lD orbitals, which is 
plotted in Fig. [SJb). Although the overall profile in Fig. 
[5Jb) is not exactly the same as that in Ref. [35|, which is 
attributed to more complicated momentum dependence 
of tunneling matrix element o 26 ' 40 ' 41 not considered here, 
it captures the increasing DOS with the magnetic field 
which is more consistent with Ref. 

The insets in Fig. [5] plot the total and quasi-lD band 
DOSs at zero field within a wider range of < 0.4. 

The peaks corresponding to the van Hove singularities 
reside at ui/ti « —0.2 and u)/t\ « 0.35, far away from the 
Fermi energy. The reason why the small energy scale like 
Zeeman energy (~ 0.003ii) can push the system to get 
closer to the van Hove singularities at energies far away 
from the Fermi energy is the help of the metamagnetism. 
In the mean-field theory, the magnetization produces an 
effective chemical potential shift as sU S" for electrons at 
orbital a and spin s. As a result, under the magnetic field 
the jump in the magnetization gives S" ~ 0.05 within 
the range of experimental interests. This leads to the 
effective chemical potential shift about ±0.18ti for U = 
3.6ti, which is large enough to push the system closer to 
the van Hove singularities. This renormalization of the 
chemical potential by the interaction is also part of the 
cause for the violation of the rigid band picture. 



D. Finite Temperature Phase Diagram 



since the changes in the Fermi surface topology inevitably 
lead to the non-monotonic behavior in pt Q t(w = 0). In 
particular, comparing Figs. Q] and 02 it is straightfor- 
ward to see that more Fermi surfaces appear near the 
(±7r, 0) and (0, ±7r) as the magnetic field increases. Since 
there are van Hove singularities near these four k points, 
Ptot(w — 0) is expected to increase. As the magnetic field 
is increased further so that the van Hove singularities are 
all covered below the Fermi surfaces, ptot{u — 0) starts 
to drop (not shown here). 

At the first glance, the entropy measurement and our 
results of the total DOS seem to contradict with the STM 
measurement. While both the entropy measurement and 
our results develop a maximum around the nematic re- 
gion in T>(ep), the STM measurement showed instead 
that T>(ep) keeps increasing even after the nematic re- 



Another intriguing observation by Rost et. al££ is the 
phase boundary of the nematic phase at the finite tem- 
perature. Since the nematic phase is bounded by the two 
first-order metamagnetic transition, the phase boundary 
can be determined by the magnetic Clausius-Clapeyron 
relation: 



AS 
AM ' 



(28) 



where /io is the permeability constant, H c and T c the crit- 
ical field and temperature, S the entropy, and M = \M\ 
the magnetization. Experimentally it was found that the 
entropy is always higher inside the nematic phase than 
the adjacent normal phases but the magnetization in- 
creases monotonically with the magnetic field. From Eq. 
[28l this implies a 'muffin'-shaped phase boundary, i.e., 
at field strengths slightly below 7.8T and above 8. IT the 
nematic phase appears at finite temperature but vanishes 
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FIG. 6: (Color online) (a) Phase diagram for higher tem- 
peratures. Magnitudes of N are represented by the color 
scales. The areas with light colors have large TV, defin- 
ing the region for the nematic order. The re-entry of the 
nematic order at higher temperature is seen at fields be- 
tween 0.0058 < naB/ti < 0.0063. (b) The DOSs of the 
all bands (solid line) and quasi- ID bands (dashed line) at 
fisB/ti — 0.006. The yellow areas refers to the energy win- 
dow bounded by -iiksT ' jt\ with temperature ksT/ti = 0.003. 
It can be seen that this thermal energy window covers a re- 
gion in which the DOS increases abruptly, driving the nematic 
phase at finite temperature. 



at zero temperature, to which we term as 're-entry' be- 
havior throughout this paper. 

By inspecting Fig. [4] closer, it can be seen that the en- 
tropy drops as entering the nematic phase from the lower- 
field boundary but raises as entering from the upper-field 
boundary. In other words, the present theory has the 
're-entry' behavior of the nematic order at the upper- 
field boundary but not at the lower-field boundary. This 
aspect is further confirmed by the temperature depen- 
dence of the nematic order parameter M as a function of 
magnetic field shown in Fig. [6[a). 

The re-entry behavior can be understood as follow- 



ing. Because the nematic phase transition is first-order, 
roughly speaking one requires U e ffT>(eF), where U e ff is 
the effective interaction for nematicity, to exceed a crit- 
ical value for the nematic ordering to occur. The mech- 
anism for nematic ordering in Si"3Rii207 based on van 
Hove singularities is all about increasing T>(ep) abruptly 
by driving the system closer to the van Hove singu- 
larities with the magnetic field. At the field strength 
slightly above the upper critical field, Viep) is large but 
still not enough for the occurrence of the nematic or- 
dering. Therefore, if the thermal energy is large enough 
to cover enough DOS within the thermal energy win- 
dow tp — ksT < €f < £f + kgT but still low enough 
so that the thermal fluctuations are small, the re-entry 
of the nematic phase at higher temperature can be pos- 
sible. As a illustration, Fig. |5[b) plots the DOS for 
/is-B/ti = 0.006 at which the nematic ordering first oc- 
curs at ksT/ti = 0.003 in our calculation. It can be 
seen that the thermal energy window for ksT/ti = 0.003 
(yellow areas) covers a region in which the DOS increases 
abruptly, consistent with the mechanism for the re-entry 
behavior discussed above. 



V. THE CASE OF THE TILTED MAGNETIC 
FIELDS (6 0) 

A. 6 dependence of the nematic ordering 

Experimentally, it has been concluded from the mag- 
netic susceptibility and resistivity measurements^^ 
that the resistive anisotropy disappears very quickly as 
the magnetic field is tilted away from the c-axis, suggest- 
ing that the nematic ordering vanishes with the increase 
of the field angle 9. Fig. presents the nematic order pa- 
rameter M as functions of [IbB and 9. It should be noted 
that since the in-plane component of the orbital Zeeman 
energy explicitly break the C4 symmetry down to the C2 
symmetry, N is non-zero as long as 9 ^ 0. Nevertheless, 
the experimentally observable nematic phase can still be 
identified by the jumps in A/". Our results showed that 
the nematicity is strongly enhanced with the increase of 
9, which does not agree with the experiments. 

The large enhancement of the nematic phase for 9 ^ 
in our calculation is due to the orbital Zeeman energy. 
The anisotropic in-plane component of the orbital Zee- 
man energy term — /is_B L x ^ a sin 9 is clearly propor- 
tional to B and largest at 6 = 90° (i.e. B || x). To 
illustrate this point, we plot the Fermi surfaces without 
any interaction for B || x with strength /xs|J5| = O.lii 
in Fig. [8l and an anisotropy can be seen. Although 
such anisotropy in the band structure is not important 
at low field, it can be amplified by the effect of the in- 
teractions, driving the system more susceptible to the 
nematic phase as the critical points is approached. As 
a result, the portion of the nematic phase in the phase 
diagram is enlarged as 9 increases from 0° to 90° as seen 
in our calculations. The disagreement between our the- 
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FIG. 7: (Color online) The nematic order parameter N as 
functions of hbB and 8. The magnitudes of Af are repre- 
sented by the color scales. For 8^0, the nematic order 
parameter is in general non-zero, but there is an area (with 
light colors) in which the nematic order parameter is sharply 
enhanced, defining the region of the nematic phase observed 
in experiments. 




FIG. 8: (Color online) The non-interacting Fermi surfaces for 
B || x and hbB = O.lti. The anisotropy in Fermi surfaces 
is already visible, especially for the areas near (±7T, 0) and 

(0,±7T). 



ory and the experimentally-determined phase diagram as 
a function of 9 could result from the strong correlation 
effect, which is beyond the scope of the mean-field theory. 



B. Anisotropy of the resistivity measurements 

While it requires more experimental efforts to confirm 
the cause of disagreement, Raghu et. al. have argued 
that the resistivity measurement might not be a good 




FIG. 9: (Color online) Illustration of the energetically-favored 
domain structures as the in-plane magnetic field is along x 
axis. When the in-plane magnetic field is weak, orbital or- 
dered phase with A/" > are dominant (shaded areas) and 
high-energy domains with Af < (green ovals) could exist, 
(b) The domain walls extend longer in y direction because 
it costs less energies if less y- than x- bonds are broken in 
a Af > background. The white oval represents the wave 
function of d yz orbital on each site. The current flows much 
more easily along y-axis than i-axis since the electrons suffer 
less domain scatterings hopping along this direction. As the 
in-plane magnetic field is strong, these high-energy domains 
vanish, thus there is no longer an easy axis for the current 
flow. 



indicator for the nematic phase in this material^ 5 -. The 
point is that the nematic phase is mostly associated with 
the states near the van Hove singularity whose Fermi ve- 
locities are too small to have sizable contribution to the 
transport properties. Therefore the observed anisotropic 
resistivity is mostly likely due to the scatterings between 
nematic domains, and the tilt of the magnetic field helps 
lie up the domains. However, if the domains are fully lied 
up, the resistivity measurement becomes insensitive to 
the nematic phase due to the diminishing scatterings be- 
tween nematic domains, despite the nematic order could 
be stronger. 

One remaining puzzle of this domain scattering argu- 
ment is why the easy axis for the current flow is per- 
pendicular to the in-plane component of the B- field 5 . 
We provide a natural argument based on the anisotropic 
spatial extension of the domain walls as explained be- 
low. Assuming the B-field lying in the xz-plane, the 
in-plane (xy) orbital Zeeman energy reads Hi n - P i ane — 
— HbB sin 9 ^2ii a Lx,iaj which couples the d xy and d xz - 
orbitals and breaks the degeneracy between the d xz and 
dj, z -orbitals. Since the d xy -orbital has lower on-site en- 
ergy due to the crystal field splitting than that of d xz , the 
<i xz -orbital bands are pushed to higher energy than d yz - 
orbital bands by this extra coupling. As a result, the ne- 
matic state with preferred dj, z -orbitals (i.e., Af > 0) has 
lower energy in the homogeneous system. At small angles 
of 9, domains with preferred d^-orbitals (i.e. Af < 0) 
could form as depicted in Fig. [9] (a) as meta-stable 
states, which occupy less volume than the majority do- 
main of Af > 0. Let us consider the shape of the domain 
walls. Because of the quasi- ID features of the d xz and 
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cL z -orbitals, the horizontal (vertical) domain wall breaks 
the bonds of the d yz (d xz )-orbital as depicted in Fig. [§] 
(b), respectively. Since the d^-orbital is preferred by 
Hi n ~piane, the horizontal domain wall costs more energy. 
Consequently, the domain structure illustrated in Fig. [9] 
(a) with longer vertical walls than the horizontal walls is 
energetically favored. 

Since the electrons suffer less domain scatterings hop- 
ping along the y-axis in this domain structure, it becomes 
the easy axis for the current flow. At large values of 
8, higher energy domains are suppressed and eventually 
vanish, and thus the resistivity measurement becomes in- 
sensitive to the nematic phase because of vanishing of the 
domain scatterings. 



C. Measurements proposed for the nematicity 

Our results have posted a possibility that the nematic 
order could occur in a larger range of the magnetic field 
for B || x than for B \\ z. Detection methods other 
than resistivity would be desirable. One feasible way is 
to measure the quasiparticle interference (QPI) in the 
spectroscopic imaging STM which has been examined in 
detail in our previous worki£. It has been predicted by us 
that if there is a nematic order, QPI spectra will manifest 
patterns breaking rotational symmetry. 

Another possible experiment is the nuclear quadruple 
resonance (NQR) measurement, which has been widely 
used to reveal ordered states in high-T c cuprates^2r— and 
recently the iron-pnictides 45 . This technique utilizes the 
feature that a nucleus with a nuclear spin I > 1 has a 
non-zero electric quadruple moment. Because the elec- 
tric quadruple moment creates energy splittings in the 
nuclear states as a electric field gradient is present, a 
phase transition could be inferred if substantial changes 
in the resonance peak are observed in the NQR measure- 
ment. Besides, since this is a local probe at the atomic 
level, it is highly sensitive to the local electronic change. 
Given that Ru atom has a nuclear spin of / = 5/ 3^ and 
the orbital ordering in the quasi-lD bands significantly 
changes the charge distribution around the nuclei, a sys- 
tematic NQR measurement as functions of magnetic field 
and field angle will reveal more conclusive information 
about nematicity. 



VI. CONCLUSIONS 

In this paper we have shown that a number of impor- 
tant properties observed in the bilayer Sr3Ru20y could 
be qualitatively consistent with a realistic tight-binding 
model together with on-site interactions treated at mean- 
field level. The band structure of this material is compli- 
cated by multibands, bilayer splitting, rotations of RuO 
octahedra, and the spin-orbit coupling, collectively lead- 
ing to the high sensitivity to the small energy scales. This 



is the main cause of the singular behavior in the evolu- 
tion of the Fermi surfaces under magnetic field. Using 
our tight-binding model and the standard mean-field ap- 
proaches on the intra- and inter- orbital Hubbard inter- 
actions, we find that for the case of magnetic field paral- 
lel to the c-axis, the nematic order, which is interpreted 
as the orbital ordering in quais-ld d yz and d xz bands, 
appears at the magnetic field at the same order of the 
experimental values ~ 8T. Moreover, we find that the 
total density of states at the Fermi energy T>(ep) under 
the magnetic field does not follow a rigid band picture, 
in agreements with the results of STM measurement^. 

The failure of following a rigid band picture is essen- 
tially a consequence of the interplay between spin-orbit 
coupling and the Zeeman energy, despite the strong cor- 
relation effect could also result in the violation of the 
rigid band shift upon doping^. Because the spin-orbit 
coupling hybridizes the quasi-lD bands and 2-D bands 
with opposite spins, the Zeeman energy naturally induces 
the reconstruction of the Fermi surfaces instead of just 
rigid chemical potential shifts. This singular behavior in 
D(ei?) also results in the divergences in the entropy and 
specific heat landscapes, since at very low temperature 
both quantities are approximately proportional to T>(ep). 
Because the divergence observed by Rost el. al^ start 
approximately at 6-7 Tesla which is not very close to the 
quantum critical point residing about 8 Tesla, a direct 
application of the quantum critical scaling seems to be 
inappropriate. The explanation of the critical exponent 
associated with this divergence could not be complete 
without taking the band structure singularity into ac- 
count in this particular material^!. 

As the magnetic field is tilted away from the c axis 
(#7^0), we find that nematic region expands instead of 
shrinking as the resistivity measurement has indicated. 
From the theoretical viewpoints, the tilt of the magnetic 
field induces an extra in-plane component of the orbital 
Zeeman energy which explicitly breaks the C4 symmetry 
down to the Ci symmetry. As argued above that this 
system is very sensitive to small energy scale, the effect 
of this extra Zeeman energy is not important at low field 
but could amplify the effect of interaction to drive the 
system toward nematicity as the quantum critical point 
is approached. As a result, the nematic phase is more fa- 
vored and stable in the presence of the in-plane magnetic 
field and it requires another Fermi surface reconstruc- 
tion at even higher magnetic field in order to weaken the 
nematic phase by reduced T>(cf). 

To explain this discrepancy between our theory and the 
resistivity measurement, we adopt the domain scattering 
argument proposed by Raghu et. aL— . Furthermore, 
we have given an explanation for another experimental 
puzzle that the easy axis for the current flow is always 
perpendicular to the in-plane magnetic field. Measure- 
ments like quasiparticle interference in the spectroscopic 
imaging STM and NQR which could detect the orbital 
ordering directly have been proposed to be more reliable 
probes for the nematicity in this material than the resis- 
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tivity measurement. 

Finally we would like to comment on limitation of the 
present theory. Although we have found the 're-entry' be- 
havior of the nematic phase, i.e., the appearance of the 
nematic phase only at the finite temperature but not at 
the zero temperature, near the upper-field boundary, the 
experiments showed this behavior near both upper- and 
lower- field boundary. In our calculations, the re-entry 
behavior is due to the increase of the density of states 
within the narrow energy window around the Fermi en- 
ergy opened by thermal energy, but we do not reject other 
schemes for the re-entry behavior. One possible scheme 
is an analogue of ferromagnetism without exchange split- 
ting proposed by Hirscb4£. He showed that the nearest 
neighbor interactions could result in a spin-dependent 
renormalization on the bandwidth (equivalently, the ef- 
fective mass). As a result, the filling for different spin 
bands can be different because of the unequal effective 
masses, leading to the ferromagnetism even without the 
exchange splitting as in the Stoner model. 

In Hirsch's original paper, the re-entry of the ferromag- 
netism at higher temperature was found. Since we only 
considered the on-site interactions in our model, such 
an effect is beyond the scope of the current theory. It 
is possible that after including the nearest neighbor in- 
teraction, the renormalizations of the bandwidths have 



novel temperature-dependences, leading to a phase dia- 
gram better consistent with the experiments. If this is 
the correct scheme, the re-entry of the nematic states 
should be accompanied by a change in the kinetic energy 
due to the effective mass renormalization, which could be 
examined by the sum rules for the optical properties^ - — . 

Another possible scheme for the re-entry behavior is 
the quantum critical fluctuations. It is well-known that 
the influences of the critical fluctuations extend from the 
quantum critical point to finite temperature in a V^-shape 
region in the phase diagram. Moreover, the critical fluc- 
tuations in this material contain not only the ferromag- 
netic but also the nematic ones. As a result, it is not 
surprising that the competition between these two types 
of critical fluctuations leads to a intriguing phase dia- 
gram at the finite temperature, and the study toward 
this direction is currently in progress. 
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